Astrophysical foregrounds and primordial tensor-to-scalar ratio constraints from 

CMB B-mode polarization observations 
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We study the efTects of astrophysical foregrounds on the abihty of Cosmic Microwave Background 
B-mode polarization experiments to constrain the primordial tensor-to-scalar ratio, r. To clean the 
foreground contributions we use parametric, maximum likelihood component separation technique, 
and consider experimental setups optimized to render a minimal level of the foreground residuals in 
the recovered CMB map. We consider nearly full-sky observations, include two diffuse foreground 
components, dust and synchrotron, and study cases with and without calibration errors, spatial 
variability of the foreground properties, and partial or complete B-mode lensing signal removal. 

In all these cases we find that in the limit of an arbitrarily low noise level and in the absence of 
the systematic effects, due to the instrument or modeling, the foreground residuals do not lead to a 
limit on the lowest detectable value of r. Nevertheless, the need to control the foreground residuals 
will play a major role in determining the minimal noise levels necessary to permit a robust detection 
of r( ^ 0.1) and therefore in optimizing and forecasting the performance of the future missions. 
For noise levels corresponding to current and proposed experiments, the foreground residuals are 
found non-negligible and potentially can significantly affect our ability to set constraints on r. We 
furthermore show how in some of these cases the constraints can be significantly improved on by 
restricting the post-component separation processing to a smaller sky area. This procedure applied 
to a case of a COrE-like satellite mission is shown to result potentially in more than an order of 
magnitude improvement in the detectable value of r. Our conclusions are found to be indepen- 
dent on the assumed overall normalization of the foregrounds and only quantitatively depend on 
specific parametrizations assumed for the foreground components. They however assume sufficient 
knowledge of the experimental bandpasses as well as foreground component scaling laws. 

PACS numbers: 98.70.Vc, 98.80.Es 



I. INTRODUCTION 

Cosmic Microwave Background (CMB) B-modc obser- 
vations are expected to be a primary source of informa- 
tion about the physics of the very early Universe, po- 
tentially providing an unambiguous proof of existence 
of the primordial gravity waves, e.g., [1-4], considered 
a telltale signature of inflation. Consequently, the CMB 
B-mode observations are a very dynamic area of the cur- 
rent research in cosmology, with multiple observatories 
being designed, built, and deployed, e.g., [5-9]. The chal- 
lenges expected in carrying out such a program are as 
impressive as its goals exciting. Besides the observation- 
specific issues, related to instrumental hardware or its 
operations, two effects: lensing-induced B-mode signal 
and astrophysical foregrounds, have come to the fore, 
both deemed as capable of setting some ultimate limi- 
tations to the exploitation of the CMB B-mode potential 
[10, 11]. In this paper we revisit these problems from a 
perspective of the parametric component separation tech- 
nique [12, 13]. Its keystone assumption is that sufficiently 
precise parametrizable frequency scaling laws are avail- 
able for each of the relevant components. At the present, 
such an assumption may look somewhat farfetched, in 
particular, in the context of the high precision required 
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for the B-mode data analysis. However, the parametric 
methods, in their modern formulation have been shown 
to perform very well at the present in the number of con- 
texts, e.g., [14-16], and this, given the theoretical and ob- 
servational work undertaken currently and aiming at un- 
derstanding all the major foreground components, seems 
to bode well for its future applications. In the context 
of this study the method provides a flexible framework, 
amenable to semianalytic, statistically robust analysis. 

In this paper we consider experimental setups opti- 
mized to ensure the lowest foreground residual level [17]. 
We use two component foreground model, including syn- 
chrotron and dust, each parametrized with one parame- 
ter, referred to as spectral indices. We note that more 
refined scaling laws with potentially more parameters can 
be straightforwardly incorporated in the formalism pre- 
sented below and would affect our conclusions only quan- 
titatively. The foreground templates used here are de- 
scribed in detail in [17, 18] and we assume nearly full 
(^ 80%) sky coverage, corresponding to the choice of 
MASK-i of [17]. Given the optimized setup and its noise 
we estimate a typical residual and compare it with the 
total statistical uncertainty. As the latter depends on the 
tensor-to-scalar ratio, r, for each value of r we determine 
respective instrumental sensitivity for which the residual 
is irrelevant given the uncertainties and compare it with 
the statistical limits due to the noise and the CMB signal 
only. 



II. METHODOLOGY 

Parametric component separation: The fiducial 
data set we consider hereafter is made of multiple-single 
frequency maps of Q and U Stokes parameters, with the 
instrumental noise assumed to be uncorrelated, both be- 
tween the pixels and channels, pixel- independent, and 
characterized by its variance, N. The corresponding 
data model we use hereafter then reads, 

dp = B {f3, ij) Sp + np = n (u) A ((3) Sp + Up, (1) 

where for each pixel p, Ais a mixing matrix parametrized 
by the spectral indices, /3, Sp - a vector of sky signals to 
be recovered and Up - instrumental noise. J7 is a pixel- 
independent, diagonal matrix with the diagonal elements, 
u;^ = i7ii corresponding to the calibration factors for each 
of the channels. The likelihood function then reads [13], 

-21n£ = ^ {dp - B Spf N-'^ {dp - B Sp) 
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where the last term is simply a prior term constraining 
the plausible values of the calibration factors. Hereafter 



we will assume that the true values of the calibration 
factors are equal to unity, Wj = 1, and that their 
uncertainty is described by an error matrix, H, which 
for simplicity is assumed to be proportional to a unit 



matrix, i.e., H^ 



(7, , ^ (5f , where cr^ is assumed not to 



depend directly on the parameters of the considered 
experiment. Moreover, throughout this paper the 
detector's bandpasses are always taken to be known 
perfectly and therefore their effects on the mixing 
matrix, B, straightforwardly calculable. The import of 
the bandpass uncertainties will be studied elsewhere. 
In the cases without calibration uncertainty, B — A 
and we simply drop the last term. In general, we will 
estimate both /3 and u: and maximize this likelihood to 
perform the component separation. 



Residual computation: The computation of the 
residuals involves two steps. First, we obtain the er- 
ror of the estimation of the spectral parameters. This 
is done using a generalization of Eq. (5) of [17], allowing 
for the calibration errors [13], and derived as the Fisher 



matrix, S, 
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of the profile likeli- 



hood, Cspec, of the likelihood in Eq. (2), i.e.. 
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which has to be evaluated at the true values of the param- 
eters, 7 = 7, where 7 stands for either /3 or u>, ,i = d/dji, 
and the matrix F defined as .F = nl^j^ ^ Sp Sp* en- 
capsulates all the information about the sky components 
needed for the parameter errors estimation. In the fol- 
lowing we will be removing the contribution to S related 
to the mode v (x [0, ..., 0, 1, ..., 1], where the zeros are as- 
signed to the spectral parameters, /3, and ones to the cal- 
ibration ones, a;, and v is normalized to one. This is done 
by replacing S by S — (f* 11v)vv*. The mode v describes 
an overall miscalibration of the final CMB map, RMS of 
which is given by a^, introducing a similar error in our 
determination of r, which typically is much smaller than 
the statistical uncertainty {Sr/r ^ 0.01 ^ crj for r ^ 0.1), 
and thus negligible. 

We use the recipe of [18] to calculate the power spectra 
of the typical noise-free foreground residuals, C^ , found 
in the separated maps, i.e.. 



k,k' j,j' 



2jkk' OLu OCu, 0<, 



(4) 



here C^ is a cross-spectrum of components i and j, and. 



OLk = 
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Residuals significance: We quantify the importance 
of the residuals as follows: 



cy-'^ 



fsky 2_^ 
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(5) 



which can be derived as a Fisher error on an overall 
amplitude, a{— 1), of a foreground template, assumed 
to be known, with the power spectrum given by C^. 
a^^ expresses statistical significance with which the 
template could be detected, had it been known, given 
the instrumental noise, Cf°^^°, and the CMB signal, 
Cf"'^{r) + 77(7)°"". Tj {< 1) denotes the fraction of 
the lensing signal left after its removal. Whenever 
a~^ is large, the residual can not be neglected in an 
analysis of the CMB map and may need to be treated 
by some additional means [19]. Otherwise, the fore- 
ground residuals will be irrelevant for the estimation of r. 

Experiment optimization: We use the approach 
described in [17] to optimize the experimental setups. 
We assume a fixed, though arbitrary, focal plane area 
during the optimization and restrict frequencies of the 
observational channel bands to range from 30 GHz to 
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FIG. 1. The significance of the foreground residuals, cr~ , Eq.(5), expected in the recovered CMB map covering ~ 80% of the sky for the cases 
with no, r] — 1.0, (left panel), partial, 77 — 0.1, (middle panel), and complete, r] — 0.0, (right panel), lensing correction, respectively. The color 
bands correspond to different calibration uncertainties as listed in the left panel with the gray color showing all the cases with (T^^j ^ after the 
removal of the mode v. The width of the shaded areas refiects the effect of varying r from 0.001, (upper edge), up to 0.1, (lower), and the dashes 
show the corresponding r — cases. The black solid lines show the case with (7^ — 0, r — 0, and rj — 1.0 as a reference. 



400 GHz. The detector noise is assumed to be constant 
in antenna temperature units. The optimization then 
tries to minimize the effective r value as proposed in [20] , 

given by J^T"" C^"™{rcB) = J2e"°" '^t- The criterion 
selection reflects the fact that we want to minimize 
the effects of the foreground residuals and thus keep 
their expected level as low as possible, irrespective of 
consequences it may have on, e.g., effective noise of the 
experimental configuration selected in such a way. The 
resulting experiment setup includes 5 frequency bands: 
V — [30,40,130,300,400] GHz occupying, respectively, a 
fraction /p(= [9, 21, 36, 25, 9] per cent) of the focal plane. 



III. DISCUSSION AND RESULTS 

Hereafter we will use the noise level of the recovered 
CMB map as a measure of the sensitivity of the consid- 
ered experimental setups. This is given by. 



"CMB 
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where we assume that CMB is the zeroth component re- 
covered in the separation procedure. The diagonal ele- 
ments of the correlation matrix, N, expressing the noise 
level of each frequency channel, can be written in antenna 
temperature units as. 
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where unet is a frequency-independent detector of 
instantaneous noise value (in /iKant-s/sec), Afp, and 
Ad{v{i)) - total and per detector effective focal plane 
area, Tobs ~ total observation time, and fip - pixel size in 
steradians. For the considered experiment we can write 



numerically, 
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The dependence of our measure of the significance of the 
foreground residuals, cr~^, on the noise level, ctcmBj is il- 
lustrated in Fig. 1, and its major features can be tracked 
back to the behavior of the parameter errors and fore- 
ground residuals as it is discussed in detail in [21]. In par- 
ticular in the low-noise regime the value of a~ ^ increases 
ex CTQjyjg whenever no calibration uncertainty is present 
or the contribution of the mode v is suppressed. This is 
due to the fact that the error on all the parameters 7 is 
driven by the first term on the right-hand side of Eq. (3), 
resulting in a self-calibrating property of the considered 
system thanks to the assumed scaling laws spanning the 
entire range of considered frequency bands. The self- 
calibration applies only to the relative calibrations fix- 
ing the calibration coefficients of the channel maps with 
precision superseding that given by the assumed priors. 
The absolute calibration of the final map is in turn al- 
ways determined by the prior term in Eq. (3) and thus 
independent on the experimental noise, as shown by the 
fiat, low-noise asymptotes of the lines, computed with the 
mode V included. For higher noise levels the calibration 
errors have significant impact on the residual level and 
should be therefore included in any meaningful analysis. 
Whenever the first term on the right-hand side of Eq. (3) 
is dominant, our results also do not depend, or depend 
only very weakly, on the foreground amplitude and on the 
observed sky area (at least as long as the foregrounds are 
nearly stationary), as the foregrounds amplitudes present 
in the expressions for S and C^ cancel. Physically, this 
means that higher levels of foreground signals lead to 
tighter constrains on their parameters, compensating for 
their higher amplitudes. 

The results from the three panels of Fig. 1 are trans- 
lated into limits on acMB, as shown in Fig. 2, by solving 




patches, for each of which we assign a different set of spec- 
tral parameters. If the patches are of roughly the same 
size, the resulting errors on the spectral parameters will 



increase approximately as 



"■pj 



leading to a tightening 
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primordial lensor-to-scalar ratio, r 

FIG. 2. Upper limits on the map noise levels, whieh ensure that the 
foreground residuals are statistieally irrelevant, are shown with solid 
lines. Each set of three lines corresponds to a different assumptions 
about the calibration errors as marked in the figure. In each set the 
lines depict the cases with no (heavy), 90% (medium), and perfect 
(thin) cleaning efficiency. The thin dashed line shows the change in the 
derived noise levels incurred as a result of restricting the sky area used 
to estimate r after the component separation step has been already 
performed. These should be compared to the thick line with ct^ — 0. 
The thick dots show the analogous noise limits based on an alternative 
criterion, r^ff: Sec. II. The shaded areas depict statistical 2a limits due 
to the noise and sky signal for three lensing cleaning efRciencies rj — 1.0, 
0.1, and 0.0 (light to dark grey). The noise levels for re-optimized 
Planck and COrE-like experiments are also shown as a reference. 



the relation, cr^^ {r, a qmb) = '^a^lcrif Hereafter, we 
use cr~^\ — 1, corresponding to a "Itr" detection of 
the residuals on the map level. In general, this value 
should be adjusted, and the curves in the figure rescaled 
by ex (Jalcrit I given a specific application envisaged for 
the output maps and 1 is used here as an illustration. For 
each r value, each curve, computed for specific assump- 
tions about the experiment and/or foregrounds, provides 
an upper limit on the experiments sensitivity so the fore- 
ground residuals will be found irrelevant for the analy- 
sis of the obtained CMB map. The gray-shaded areas 
show the statistical uncertainties, corresponding to a dif- 
ferent level of gravitational lensing signal cleaning. We 
note that the foreground residual limits do not prevent 
detecting arbitrarily low value of r assuming that a suffi- 
ciently sensitive observation can be performed. Instead, 
the lower limit on r can arise due to a residual level of the 
lensing-induced B-mode signal left over from some clean- 
ing procedure [22-24]. This remains true when the cali- 
bration errors are included but also when the spatial vari- 
ability of the foregrounds is allowed for, and will hold at 
least as long as no significant deviation from the assumed 
component scaling laws is observed. To see the effects 
of the spatial variability of the spectral indices we as- 
sume that the sky is subdivided into rip non-overlapping 



of the noise constraints in Fig. 2 by the same factor. For 
comparison the (orange) line in Fig. 2, labeled "spatially 
varying Pdust^ , shows a result of implementing the Stol- 
yarov approach [25], which also leads to more restrictive 
noise constraints, but without introducing an ultimate 
limit on r. 

We also note that by decreasing the statistical un- 
certainty of the map we increase cr~^, as the residual 
becomes easier to be spotted, and thus the requirements 
on the noise need to be tighter to ensure that the fore- 
ground level is decreased accordingly. This, for instance, 
explains why any lensing cleaning in Fig. 2 renders a 
tighter limit on the noise. Conversely, re-sorting for 
the r estimation to a smaller map of the sky, than 
what has been used for the component separation, will 
increase the variance and lower a~^ , allowing us to tune 
appropriately the sky area to extend the range of de- 
tectable values of r given a fixed instrumental sensitivity. 
This will result in lower statistical significance of the 
detection but will ensure that bias is negligible. This is 
illustrated in Fig. 2 where the (magenta) arrows show a 
change in the noise upper limit in the perfect calibration 
case, a^j = 0, with no lensing cleaning, rj = 1.0, due to 
using on the second step only half of the area of 80% 
of the full sky as used for the component separation. 
This, for a COrE-like experiment (i.e., the proposed 
COrE experiment [7] optimized as described earlier), 
could extend its capability to detect r reliably down to 
8 X 10~^ (20"), what could be compared to r ^ 4 x 10^"^ 
limit (2cr) potentially achievable, if the foregrounds were 
absent. We note that the trimming can be made even 
more efficient if the retained sky is selected to ensure 
the lowest possible foreground amplitude. If no extra 
trimming is done, then given our criterion for a^^ the 
COrE-like lower limit on r is found to be r ^ 3 x 10^^, 
what is at least formally within reach of a suborbital ob- 
servation with similar sensitivity per pixel but observing 
0(1)% of the sky [18, 19]. The statistical significance 
of the former limit is '^ 25cr, (vs. 2(t in the suborbital 
case) indicating that the experimental sensitivity of such 
observations should be driven by the foreground separa- 
tion, not by statistical uncertainties only, but also that 
a further improvement of the limit on r could be plausi- 
ble if extra assumptions and processing are included [19]. 

The results obtained here demonstrate that in an ab- 
sence of such post-component separation processing and 
with calibration uncertainties as typically present in ac- 
tual experiments the noise levels required for an unam- 
biguous and robust determination of r are on order of 
C'(10~"'^)/zKarcmin, significantly below the noise levels 
for the currently considered satellite mission concepts. 
Moreover, if the lensing contribution left over after its 
cleaning is higher than ^ 10% of its initial value, the 



dependence of the noise levels on the targeted value of 
r is rather weak. This emphasizes that once the suffi- 
cient noise level is indeed attained the measurable values 
of r would be limited only by the statistical uncertain- 
ties. On the contrary, a failure to reach such a noise level 
may render the experiment incapable of setting any con- 
straints on r of current interest. If the lensing could be 
cleaned nearly perfectly, rj ;S 10%, lower noise levels lead 
to a progressively lower limit on the detectable r. 

Summarizing, we have studied the importance of the 
foreground residuals left over from the maximum likeli- 
hood parametric component separation procedure on the 



detection of the primordial tensor-to-scalar ratio coeffi- 
cient, r, by nearly full-sky CMB B-mode experiments. 
We have found that though the foreground residuals are 
likely to be a major driver in defining the sensitivity re- 
quirements for such experiments, they do not on their 
own lead to any fundamental lower limits on detectable 
r, at least as long as sufficiently precise frequency scal- 
ing models are available. These will be rather set by the 
uncertainty due to the lensing signal present in the maps 
after its cleaning. We note that the latter may also in 
turn depend on the presence of foregrounds and instru- 
mental noise [26, 27], an issue we address elsewhere [21]. 



[1] U. Seljak and M. Zaldarriaga, Physical Review Letters 

78, 2054 (1997), arXiv:astro-ph/9609169. [15 

[2] W. Hu and M. White, New A 2, 323 (1997), arXiv:astro- 
ph/9706147. 

[3] M. Kamionkowski, A. Kosowsky, and A. Stebbins, 

Physical Review Letters 78, 2058 (1997), arXiv:astro- [16 
ph/9609132. 

[4] M. Kamionkowski and A. Kosowsky, Phys. Rev. D 57, 

685 (1998), arXiv:astro-ph/9705219. [17 

[5] K. Arnold, P. A. R. Ada, A. E. Anthony, F. Aubin, 

D. Boettger, J. Borrill, C. Cantalupo, M. A. Dobbs, [18 
J. Errard, D. Flanigan, et al., in Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference 
Series (2010), vol. 7741 of Society of Photo- Optical In- [19 
strumentation Engineers (SPIE) Conference Series. 

[6] B. Reichborn-Kjennerud, A. M. Aboobaker, P. Ade, 

F. Aubin, C. Baccigalupi, C. Bao, J. Borrill, C. Can- [20 
talupo, D. Chapman, J. Didier, et al., in Society of 
Photo-Optical Instrumentation Engineers (SPIE) Con- 
ference Series (2010), vol. 7741 of Society of Photo- [21 
Optical Instrumentation Engineers (SPIE) Conference [22 
Series, 1007.3672. 

[7] The COrE Collaboration, ArXiv e-prints (2011), [23 
1102.2181. 

[8] J. Aguirre, A. Amblard, A. Ashoorioon, C. Baccigalupi, 

A. Balbi, J. Bartlett, N. Bartolo, D. Benford, M. Birkin- [24 
shaw, J. Bock, et al., ArXiv e-prints (2009), 0903.0902. 

[9] A. Kogut, D. J. Fixsen, D. T. Chuss, J. Dotson, E. Dwek, [25 
M. Halpern, G. F. Hinshaw, S. M. Meyer, S. H. Moseley, 
M. D. Seiffert, et al., ArXiv e-prints (2011), 1105.2044. 
[10] M. Tucci, E. Martinez-Gonzalez, P. Vielva, and J. De- [26 
labrouille, MNRAS 360, 935 (2005), arXiv:astro- 
ph/0411567. [27 

[11] M. Betoule, E. Pierpaoli, J. Delabrouille, M. Le Jeune, 

and J. Cardoso, A&A 503, 691 (2009), 0901.1056. 
[12] W. N. Brandt, C. R. Lawrence, A. C. S. Readhead, J. N. 

Pakianathan, and T. M. Fiola, ApJ 424, 1 (1994). 
[13] R. Stompor, S. Leach, F. Stivoli, and C. Baccigalupi, 

MNRAS 392, 216 (2009), 0804.2645. 
[14] H. K. Erikson, A. J. Banday, K. M. Gorski, and P. B. 



Lilje, ApJ 612, 633 (2004), arXiv:astro-ph/0403098. 
H. K. Eriksen, C. Dickinson, C. R. Lawrence, C. Bacci- 
galupi, A. J. Banday, K. M. Gorski, F. K. Hansen, P. B. 
Lilje, E. Pierpaoh, M. D. Seiffert, et al., ApJ 641, 665 
(2006), arXiv:astro-ph/0508268. 

J. Dunkley, E. Komatsu, M. R. Nolta, D. N. Spergel, 
D. Larson, G. Hinshaw, L. Page, C. L. Bennett, B. Gold, 
N. Jarosik, et al., ApJS 180, 306 (2009), 0803.0586. 
J. Errard, F. Stivoh, and R. Stompor, Phys. Rev. D 84, 
063005 (2011). 

F. Stivoh, J. Grain, S. M. Leach, M. Tristram, C. Bac- 
cigalupi, and R. Stompor, MNRAS 408, 2319 (2010), 
1004.4756. 

Y. Fantaye, F. Stivoli, J. Grain, S. Leach, M. Tristram, 
C. Baccigalupi, and R. Stompor, Journal of Cosmology 
and Astroparticle Physics 2011, 001 (2011). 
A. Amblard, A. Cooray, and M. Kaplinghat, 
Phys. Rev. D 75, 083508 (2007), arXiv:astro- 
ph/0610829. 

J. Errard, G. Fabbian, and R. Stompor (in preparation). 
L. Knox and Y.-S. Song, Physical Review Letters 89, 
011303 (2002), arXiv:astro-ph/0202286. 
M. Kesden, A. Cooray, and M. Kamionkowski, Phys- 
ical Review Letters 89, 011304 (2002), arXiv:astro- 
ph/0202434. 

U. Seljak and C. M. Hirata, Phys. Rev. D 69, 043005 
(2004), arXiv:astro-ph/0310163. 

V. Stolyarov, M. P. Hobson, A. N. Lasenby, and 
R. B. Barreiro, MNRAS 357, 145 (2005), arXiv:astro- 
ph/0405494. 

W. Hu and T. Okamoto, ApJ 574, 566 (2002), 
arXiv:astro-ph/0111606. 

K. M. Smith, A. Cooray, S. Das, O. Dore, D. Han- 
son, C. Hirata, M. Kaplinghat, B. Keating, M. Loverde, 
N. Miller, et al., in American Institute of Physics Con- 
ference Series, edited by S. Dodelson, D. Baumann, 
A. Cooray, J. Dunkley, A. Fraisse, M. G. Jackson, 
A. Kogut, L. Krauss, M. Zaldarriaga, & K. Smith (2009), 
vol. 1141 of ylmerican Institute of Physics Conference Se- 
ries, pp. 121-178, 0811.3916. 



